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Abstract 

In nature, barchan dunes typically exist as members of larger fields that display enig- 
matic structures that cannot be readily explained by external forcing. To explore the 
possibility that observed patterns self-organize, we built a numerical model that treats 
barchans as discrete entities that obey known empirical relationships and interact with 
one another according to simplified rules derived from theoretical and numerical work, 
and from field observations. A rich array of patterns, similar to those observed in na- 
ture, emerge from these relatively simple behaviors and interactions, offering a cohesive 
depiction and a potential explanation of field-scale phenomena. This is the first model to 
integrate calving and results support the hypothesis that it exerts a first order control on 
the system; It is fundamental to the formation and stability of field patterns and when 
the process is disabled, all dunes in all runs grow without bound on short time-scales. 
Insights derived from this synthesis of dune-scale phenomena provide testable predictions 
for future empirical work. 

1 Introduction 

Barchans are mobile, crescent-shaped dunes that form ubiquitously atop hard, fiat surfaces 
wherever an approximately unidirectional fluid flow transports a limited amount of bedload 
sediment. At the single dune-scale, coupled models of sediment transport and fluid dynamics 
explain their morphodynamic origination and reproduce their characteristic morphology and 
kinematics: Barchans emerge from the instability of a flat sand bed and march incrementally 
forward (via avalanching) in the direction of the mean flow at a rate inversely related to 
their size (Kroy et al., 2002; Duran et al., 2010). Barchans change size as a function of 
the relative rate that sand is obtained across their upwind-facing 'back' and lost from their 
downwind-facing 'horns'. Given the way horn and dune width scale, larger barchans lose 
proportionally less sand. As a result, barchans possess a morphometric instability - dunes 
above a certain size should tend to grow without bound, while those below should tend to 
shrink and disappear - that is somehow overcome in nature, as barchans appear to persist 
while maintaining approximately constant sizes (Hersen et al., 2004; Duran et al., 2010). 

Barchans also rarely exist in isolation but decorate desert floors, continental shelves, 
neves, and celestial bodies en masse, as members of larger fields. Despite differing granular 
materials and fluid mediums, barchans from these diverse environments share more than their 
mere characteristic crescent shape: The fields themselves appear strikingly similar. Fields 
are typically heterogeneous, featuring large areas where dunes are concentrated and others 
where they are absent. Where concentrated, barchans sometimes cluster into distinguishable 
patches wherein dunes align themselves horn-to-back, resembling 'birds fiying' in a convoy. 



Other times barchans congregate in expansive downwind 'corridors' where dunes maintain a 
roughly similar size and spacing. However, this approximate homogeneity can be interrupted 
by sections of bare ground or dunes of outlier size. In addition, extensive dune fields can 
themselves be composed of individual corridors with different characteristic dune sizes and 
spacing. This enigmatic structuring, with heterogeneities and homogeneities on different 
scales and of diverse types, is not readily explained by external forcing (Elberhiti et al., 
2008). 

Forays into dune-scale dynamics, the time-scale commensurate with human observation, 
have contributed insights pertinent to the grander geologic evolution of barchan fields. How 
these smaller-scale, well-studied phenomena culminate into the larger-scale patterns charac- 
teristic of barchans fields, however, remains unresolved. Internal dynamics are thought to be 
responsible (Lima et al., 2002; Kocurek et al., 2010; Genois et al., 2012) and the hypothesis 
that calving, whereby larger dunes eject smaller dunes, is fundamental has also been posited 
(Elbelrhiti et al., 2005) but not yet tested. To explore the self-organization of field-scale 
patterns (Baas, 2007), we built a numerical model that treats barchans as discrete entities 
that obey known empirical relationships and interact according to simplified rules derived 
from theoretical and numerical work, and from field observations. Parameterizations herein 
(for numeric values not reported in main text, see Supplementary Material) are tailored to 
the sub-aerial barchans near Tarfaya in the Western Sahara based on extensive previous work 
at the dune-scale (Andreotti et al., 2002a, 2002b) and on field-scale structuring (Elbelrhiti 
et al., 2008). Our goal is to present potential insights about field-scale organization, so we 
qualitatively describe our findings and save statistical characterizations of output for ensuing 
efforts. 

2 Model 

Barchans are typically separated from one another by bare ground so we model them as 
individual objects, each with a finite cross-wind width, W, down-wind length, L, and volume 
approximated as a pyramid, V = W^/AO (Elbelrhiti et al., 2008). In plan view, we partition 
W into the characteristic features of barchans morphology: A central body defined by a slip 
face width, Ws, that is flanked by horns with a combined width, Wh, 

Wh = A + aW, (1) 

where A and a are numerically derived constants (Hersen et al., 2004). An initially empty 
field with a crosswind width, Y » W, and down-wind length, X W, is resolved into square 
cells with sides c = Wo/4, where Wq, the width of an elemental dune, is equal to the fastest 
growing wavelength that arises from the instability of a nearly flat sand bed. Sediment flux is 
partitioned into two different modes: Free flux, qf, represents wind-blown sand, and bulk flux, 
Qb, represents sand transported as migrating dunes. The algorithm below is iterated every 
time step, dt, where dt ^ T, where T is a characteristic time related to how long it takes 
to transport the volume of sand equivalent to an elementary-dune volume, Vq: T = Wg/gsat; 
where (/sat is the 'saturated flux' that a given wind regime will transport across a flat, sand 
covered surface (Andreotti, 2002a). (For sub-aerial dunes in the Tarfaya field, Wq ~ 20 m 
and T c± 2000 days, see Supplementary Material). 



2 



2.1 Sand Supplied 



At the upwind boundary, the field is suppHed with a total sand flux, qt^i, consisting of a 
'free flux' qf^i (uniformly distributed along Y), and qb^i, a 'bulk flux' which represents the 
spontaneous formation of uq elemental dunes (Elbelrhiti, 2012), 

no = ^ dt, (2) 

positioned randomly along Y. 



2.2 Sand Transported 

All cells are then assigned a free flux, qf. Traveling considerably faster than dunes, we prop- 
agate Qf^i instantaneously and directly downwind (neglecting lateral diffusion for simplicity) 
so qf = qf^i until a dune is encountered, and then is modifled to reflect how dune morphology 
impacts wind-blown transport: Downwind of a slipface, = to mimic sand trapping due 
to boundary layer separation, while downwind of each horn, qj = Qsat to capture sand leaking 
at the saturated rate (Hersen et al., 2004). When a cell is occupied by only part of a horn or 
a slip- face, (//is calculated as a weighted average, dictated by the portion of the cell occupied 
by either feature. For bulk transport, we assume no lateral migration and advance all dunes 
a downwind distance inversely related to W through the Bagnold-like translation speed, v, 

where 6 is a numeric constant and Ac is the cut off wavelength below which dunes disappear 
(Elbelrhiti et al., 2008 and references therein). 



2.3 Dunes Change Volume 

The net sand gain/loss of a dune is calculated as the difference in sand received at a rate of 
qf across W, and lost at the saturated rate across both horns: 

^ = qfW - qsa.tWh- (4) 

(When dunes shrink according to Q to the width Ac, they are assumed to lose height until 
the slip face disappears; dunes below Ac in the model lose sand at the saturated- flux rate 
across their whole width.) Since barchans introduce spatiotemporal heterogeneities in g/, the 
location of a dune relative to upwind dunes has bearing on the flrst term in Q. 



2.4 Dunes Collide 

We allow for two possible outcomes when dunes collide (Hersen and Douady 2005; Katsuki 
et al., 2005). If the center of mass of a dune falls within the body of another dune the two 
colliding dunes coalesce into a single composite dune that is located at their volume weighted 
average position. If two dunes touch but the above condition is not met, for example contact 
occurs only between their horns, they remain separate. 
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2.5 Dunes Calve 

The surface of a dune is subject to the same instabihty affecting a flat sand bed (Elberhibiti 
et al., 2005), so ah modeled dunes sport superimposed bedforms. Based on observations, 
surface perturbations with initial heights, /iq, form with a downwind spacing of Wq atop both 
dune flanks, propagate with a speed calculated from ([s]), grow exponentially at a rate, a, and 
are 'calved' as elemental dunes if they develop a slip face (sufficient aspect ratio, h{t)/Wo) 
by the time they reach the end of their host horn (Elbelrhiti et al., 2005). This formulation 
allows only dunes above a certain size, Wc, to produce calves. Yet dunes in higher density 
areas appear to begin calving at smaller sizes: The turbulent wake of one dune exaggerates 
perturbations on a nearby downstream dune. We account for this by letting Hq assume one of 
two different values depending on whether a dune is located inside or outside the wake zone 
of an upwind neighbor (Baddock et al., 2007). We subsequently refer to the two different 
values of Wc as the 'un-crowded' and 'crowded' calving widths. All successfully calved dunes 
are located immediately downwind from their parent horn and are treated akin to all other 
dunes in subsequent time steps. 

For dunes with widths well above Wc, calving frequency depends on how quickly successive 
waves propagate to the end of their host horn. However, as we describe in the Results and 
Discussion section, calving frequency within a field is effectively limited by the rate of sediment 
input into a dune; a calving event reduces W < Wc, requiring a period of subsequent dune 
growth before the next calving event can occur. 



3 Results and discussion 

Initial experiments revealed that output varies dramatically with total sand supply (Ewing 
and Kocurek, 2010; Eastwood et al., 2011) but is rather insensitive to the partitioning between 
the two transport modes (see Supplementary Material). To clearly distinguish between the 
consequences of external forcing and internal dynamics we hold qf^i at zero and vary only qt,^i, 
meaning we change the rate of dune seeding between simulations. 

A diverse array of patterns, similar to those observed in nature, result as model output 
(Figure [T]), consistent with the conception of barchan fields as emergent phenomena (Baas, 
2007): Simulated dunes appear in distinct spatial arrangements as a result of small-scale 
processes in which the large-scale effects are not obvious. The negative sand balance imposed 
on all dunes seeded at the upwind boundary, across all runs, forces any individual in isolation 
to vanish. Therefore, interactions between dunes are crucial to the formation of any non- 
empty model dune field. Furthermore, this imbues stochasticity to the sand supply: By 
releasing sand as they vanish, seeded dunes function as spatially concentrated sources of qf 
that can feed any dune, 'luckily' located directly downwind. Setting dV/dt = in (|4| reveals 
one unstable equilibrium dune width. Woo, 

Woo= ^--^^ , (5) 

that depends on qj (Hersen et al., 2004). This spatial variation in qf in turn leads to more 
than one Woo, potentially permitting divergent behavior of equivalently sized dunes; 'lucky' 
ones, with high qf from upwind, can grow, while others shrink and disappear. 

Internal dynamics, while fundamental, appear inextricably coupled to external forcing: 
The sand supply influences individual dune behavior, with higher input fluxes leading to a 
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Figure 1: Visual appearance of steady state fields varies with the sand supply rate at the up- 
wind boundary (A) q^^i = 0.20 gsat, (B) qh,i = 0.25 gsat, (C) qb,i = 0.30 gsat- As parameterized 
for the Tarfaya Field, where Wq = 20 m and T = 2000 days, these results represent a 5 km by 
30 km field after ~ 5000 years. (For time series depicting approach and maintaince of steady 
state, see Supplementary Material.) Model structures compared natural dune configurations 
observed in the Tarfaya Field from Google Earth: (A) Dune Lines, (B) 'Ducks-Flying', (C) 
Intra-dune hollow area, (D) Dune diagnols. 

higher proportion of surviving seeded dunes, for example. Individual dune behavior, in turn, 
determines how the interaction processes unfold to ultimately control field-scale structuring. 
Results highlight the overarching importance of dune calving, so we first discuss its role in 
the formation of specific patterns, and then elaborate these findings in a section dedicated to 
how the process plays out in the context of a field. 

3.1 Sand Supply 

There is a threshold sand supply, q^^i, below which simulated fields remain vacant (which varies 
given the probabilistic behavior of this model, as larger domains and/or longer durations both 
make low-frequency occurrences more likely to be observed) and beyond which fields become 
increasingly more populated (Figure [T]). 

3.2 Low Sand-Flux and 'Birds-Flying' Pattern 

When gfe.j is relatively low, fields are sparse and dunes appear in spatially distinct patches 
(Figure [i]A) . These emergent structures resemble 'birds- flying', exhibiting a V-shape centered 
on one, upwind dune and with all other patch members aligned horn-to-back (Figure [ip). 
As qh^i increases between runs and as qf increases downwind within a given run (due to 
disappearing dunes), patches become more prevalent and individual patches also tend to host 
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more dunes. This observed non-linear response to sand supply at the field-scale follows from 
single dune-scale dynamics as well as from how dunes operate as a collective, within a patch. 

The single dune located at the upwind apex of a patch initiated its formation and subse- 
quently regulates its expansion. Patches begin when a dune stochastically captures enough 
sand from disappearing upwind dunes to surmount Woo- This pioneer then grows to the 
calving width, Wc, and ejects calves. However, the sand balance of these progeny is negative 
(see Supplementary Material), implying that calf survival is tied to stochastic fluctuations in 
sand flux in a manner similar to seeded dunes. 

Yet unlike a seeded dune, a calf is mostly sheltered from variations in qj by its parent. 
Instead of buffering against fluctuations, however, a parent transmits them downwind by 
calving intermittently. Once a parent calves again, this subsequent calf disrupts the sand 
supply of the extant calf, increasing its net rate of sand loss (see Supplementary Material). 
While an upwind calf accentuates the sand-deficit of a down-wind calf, their relative spacing 
and sizes sometimes allows for their coalescence. This enables some calves to grow downwind 
of a pioneer horn, slowing their migration and allowing them to capture additional calves. If 
they become sufficiently sized, they too begin calving and some of their progeny, likewise, go 
on to form another generation centered on their horns. In this manner, model patches expand 
laterally and elongate in a stochastic fashion as they travel downwind. 

The rate of patch initiation and expansion is controlled by qb^i- As qt^i increases, the 
probability that seeded dunes are successfully established increases. Conceptually, this is a 
conditional probability where the rate of dune injection influences the probability, P, that an 
isolated dune gains an upwind neighbor, 

W 

The additional dependence of P on reveals a positive feedback: If a dune acquires an 
upwind neighbor, it grows and increases its chances of acquiring another upwind neighbor and 
growing even further. As a consequence, as q\,^i increases more dunes successfully surmount 
Woo and proceed to found more rapidly expanding patches. 



3.3 High Flux Runs and Corridor-like Structuring 

As q^^i increases further, fields become more homogeneous in the sense that dunes are observed 
throughout the domain (Figure [TJB,C), and locally within a given region, display a charac- 
teristic size (see Supplementary Material). Simulated dunes also appear in distinct spatial 
arrangements observed in natural fields, for example, in lines and diagonals (Figure [lJi],F). 
Dune density is sufficient to essentially guarantee that any given dune has an upwind neighbor 
and is therefore interacting with another dune. Bare areas intersperse the otherwise contin- 
uous sprinkling of dunes so these fields also possesses some degree of spatial heterogeneity 
(Figure [Tp). As seen in lower flux runs, this emerges despite the uniform sand input so this 
result is consistent with other research reporting that intra-field hollows and occupied areas 
exhibit no notable differences in external forcing (Elbelrhiti et al., 2008). Given this display 
of features also common to natural barchan corridors, we refer to the patterns resulting from 
these higher-flux runs as simulated corridors. 

On occasion, a dune grows large enough to substantially slow its migration and instigate 
a positive feedback: It then cannibalizes upwind dunes and interrupts the otherwise homo- 
geneous pattern in its vicinity by starving downwind dunes of sand and creating a vacant. 
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Figure 2: The dependence of field patterns on the sand supply rate, (A) g^^j = 0.25 Qsat and 
(B) Qb^i = 0.40 gsat, diminshes with distance downwind: Fields appear to approach a common 
attractor. (For histograms demonstrating local dune sizes in regions 1-4, see Supplementary 
Material.) As parameterized for the Tarfalya corridor, these results represent a 5 km by 
100 km field after ~ 10000 years (~ 1800T). 



downwind area (see Supplementary Material). This infrequently realized event is exaggerated 
by our conservative formulation of the calving process. In nature, mega-dunes are sufficiently 
wide that multiple wave trains flank each horn (Elbelrhiti et al., 2008). This presumably 
increases their calving rate and if accounted for would retard such runaway growth. Further- 
more, adding a size constraint imposed by the Atmospheric Boundary Layer (Andreotti et 
al., 2009) would also place an upper limit on the size of these dunes. 

3.4 Upwind Dependence on Input Flux, Downwind Independence 

After an initial spin up time (determined by how long it takes dunes to migrate to the 
downwind end of the domain) all simulated fields attain a statistically steady state (with 
the exception of the occasional development of mega-dunes, see Supplementary Material). 
Despite persisting differences in the upwind portion of the domain for runs with different qb^i, 
downwind portions of fields are strikingly similar (Figure [2]): The dune-dune interactions, 
when spatially iterated between successive dunes and collections of dunes over a sufficient 
distance, appear to lead toward a common attractor. Although most natural fields lack 
these expansive downwind lengths, convergent behavior could nonetheless be tested for with 
empirical work focused on down- wind field trends which, unlike transverse structure (Elberhiti 
et al., 2008), have not yet been systematically examined. 

3.5 Dune Calving Within a Field 

This work supports the hypothesis that calving plays a fundamental role in field-scale phe- 
nomena (Elbelrhiti et al., 2005). When the process is disabled, dunes grow without bound 
and fields lack any resemblance to natural ones (see Supplementary Material). When en- 
abled, calving prevents this runaway growth and allows patterns to develop which maintain 
a roughly realistic appearance. 

Model dunes appear to overcome their morphometric instability predicated by Q by 
calving, a sand loss term unaccounted for by the equation. As a growing dune reaches Wc 
and ejects a calf, its volume is diminished and its width drops below Wc, interrupting the 
process: The next wave to reach the horn tip will not quite develop a slip-face and therefore is 
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not ejected. However, the sand loss associated with one calf is not significant enough to drop 
its size below Woo (given our parameterizations, Woo Wc, see Supplementary Material). 
After recovering enough sand to again surmount Wc, it then calves and ceases calving once 
again. In the model, this negative feedback keeps dunes toggled between a narrow width range 
and leads to what, locally at the community scale, is a 'single' stable size that is approximately 

Wc. 

Moreover, the time lapsed between two successive calving events depends on W: Larger 
dunes recover sand more quickly ^ and therefore also release it more rapidly downwind 
by calving at a higher rate. Additionally, since the sand supply within a field is inherently 
stochastic, calving is also stochastic and is moderated by the behavior of neighboring dunes. 
For example, noting that dunes are typically centered on the horn of an upwind neighbor, 
a dune can absorb a calf too - which would cause a dune just below Wc to become large 
enough to calf again. Therefore Q is only loosely applicable to dunes in structured fields 
where the rate of sand gains and losses appears coupled. Short-term empirical measurements 
grounded our formulation of the calving process (Elbelrhiti et al., 2005) and longer-term 
empirical observations, for example on the relationship between calving frequency and dune 
size and field density, is needed to evaluate whether or not this calving-frequency-regulation 
mechanism occurs in natural dune fields too. 

The correspondence between the characteristic dune size of a region and Wc - either the 
crowded or un-crowded value (Figure [2| see Supplementary Material) - also points to calving 
as the model mechanism regulating dune size. Moreover, the emergent, inverse relationship 
between the size and density is consistent with the structuring observed in the Tarfaya dune 
field (Elbelrhiti et al., 2008). If this mechanism regulates dune size in nature, it is still unclear 
how different corridors hosting different sized dunes, and therefore by inference different 
calving sizes, can be juxtaposed within the same field. The only plausible explanation offered 
by this model is regional variation in sand supply: Areas with higher regional sand supplies 
could trigger the formation of a corridor of densely spaced dunes with sizes clustering around 
the crowded calving width, while areas with lower local supplies would feature dunes closer 
to the un-crowded width or organized into patches. This suggestion is consistent with field 
data (Elbelrhiti et al., 2008) if it is reinterpreted in light of the non-linear response of the 
model system to variations in input flux: The differences in regional sand fluxes within the 
Tarfaya field (~ O.OS^sat) ai"e significant enough to cause different characteristic dune sizes in 
the model, near the sand source (i.e. in the upwind portions of the modeled fields) and could 
result from regional variations in lithology, topography, or fluxes at the field boundaries. 

3.6 Sensitivity to Parameters and Robustness of Potential Insights 

Model results depend in detail on value of several parameters and process representations. 
Although some of the parameterizations in this initial modeling effort are not tightly con- 
strained, the results presented here arise from what we believe are the most reasonable choices 
among the simple formulations we have experimented with. In addition, as long as calving 
is incorporated in the model, varying most parameters alters the range of input fluxes that 
produce the various patterns and behaviors, but does not fundamentally change the qualita- 
tive results we have emphasized here. In particular, the emergence from internal dynamics of 
spatial heterogeneity within a field - including patchiness at relatively low input fiuxes and 
variations in dune density with higher sediment supplies - results robustly from combinations 
of parameterizations we have tried. As a next step, establishing and improving the quantita- 
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tive match between model and natural patterns would require integrating empirically based 
and calibrated formulations for model processes, especially dune collisions and calving, which 
here are highly simplified in this initial work. 



We thank ANR Zephyr grant #ERCS07_18 for funding. 
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A Supplementary Material 



A.l Model 

Numeric values of parameterizations used to simulate Western Sahara barchans, listed in 
order of appearance in text and drawn from Hersen et al. (2004), Andreotti et al. (2002a, 
2002b), Elbelrhiti et al. (2005), and field observations of B. Andreotti (denoted by *) 



Parameter 


Numeric value 


a 


0.05 


A 


4.65 m 


Wo 


20 m 


QsELt 


0.2 mVday 


Woo 


~ 30 m ~ 1.5Wo 


(from assuming qf = 0.2gsat) 






2000 days 


dt, where dt <^T 


1 day 


b 


45 


Ac 


16.6 m ~ 0.75Wo 


a 


0.015 m/day 


ho, outside wake zone* 


0.05 m, so Wc ~ 100 m ~ 5Wo 


ho, inside wake zone* 


0.165 m, so Wc ~ 60 m ~ SWq 


h{t)/Wo (aspect ratio) 


0.1 m ~ O.OSWo 



A. 2 Results and discussion 

Analytically, calved dunes evolve volumetrically according to, 

dV 

— = 0.5fet[A + aWc] - gsat[A + aWo], (7) 

where the first term in Q has been modified to represent sand supplied to a calf by one horn 
of its parent. Setting dV/dt = and solving for the parent width that produces stable calves 
(Wc 130 m ~ 6.5VFo) implies that parent dunes, as parameterized, cannot support a calf 
through sand-flux exchange alone. 

The first term in ([T]) represents the amount of sand leaked to a calf from a parent. After 
that parent has a subsequent calf, the amount of sand leaked downwind to its earlier calf, is 
simply the rate at which the calf loses sand, or gsat[A + aWo]i which is less than the first term 
and leads to the conclusion that subsequent calves detrimentally disrupt the sand supply of 
an earlier calf. 

In regards to the stochastic, run-away growth of large dunes (Figures [6]|8j), it should be 



noted that in nature, unusually large 'mega dunes' can occur within a dune field. These large 
barchans are sometimes assumed to form at the upwind boundary and then are used to back 
calculate the relative age of the dune field. This estimation then provides a rough proxy for 
the length of time that external conditions have been constant (Hesse, 2009). However, in 
the model these large dunes form rapidly in place, far downwind of the upwind boundary, 
suggesting that large-dune migration rates and positions could be misleading indicators of 
field age. If field-age estimates are systematically too large, this result may have implications 
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for paleo-environmental interpretation, especially for extra-terrestrial settings where other 
proxies are unavailable. 

Additional Reference (cited here in Supplementary Material but not in main text) 
Hesse R (2009) Do swarms of migrating barchans dunes record palcocnvironmcntal changes? 
- A case study spanning the middle to late Holocene in the Pampa de Jaguay, southern Peru. 
Geomorphology 104, 185-190. 
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Figure 3: Demonstration of the insensitivity of model output to flux partitioning input for 
lower flux runs where qt^i = 0.20gsat- Note that the qualitative patchy pattern is unchanged 
where (A) qj^i = O.lOgsat and Qb^i = O.lOgsat (B), qf^i = 0.05gsat and qb^i = O.lSgsat and (C) 
qf^i = O.Ogsat and qb,i = 0.20q'sat- 
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Figure 4: Demonstration of the insensitivity of model output to flux partitioning input for 
higher flux runs where qt^i = 0.30gsat- Note that the qualitative corridor pattern is unchanged 
where (A) qf^i = O.lOgsat and QbA = 0.20gsat (B), qf^i = 0.05gsat and qt^ = 0.25gsat and (C) 
qj^i = O.Ogsat and qb,i = O.SOfet- 
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Figure 5: Time series revealing model approach (A) and maintenance (B,C) of a statistical 
steady state for a lower-flux run, where qhi = 0.20(;sat- 
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Figure 6: Time series revealing model approach (A) and maintenance (B,C) of a statistical 
steady state for a higher-flux run, where qb^i = 0.25gsat 
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Figure 7: Example of stochastic and runaway growth of a mega-dune in a lower-flux run, 
where qb^i = 0.20gsat- Note this infrequent event happens far from the upwind boundary. 
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Figure 8: Example of stochastic and runaway growth of a mega-dune in a higher-flux run, 
where qb^i = 0.40gsat- 



18 



t 



i 



A. t = loor 



c 



B. t=i35r 



125Wo 



500Wo 



It t 



_s c £ E 



C. t= nor 

E 



1250Wo 



1250Wo 



Figure 9: Disabling calving process leading to run away growth of all dunes on short time 
scales, example from lower-flux run, q\)^i = 0.20gsat- 
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Figure 10: Disabling calving process leading to run away growth of all dunes on short time 
scales, example from higher-flux run, qt^j = 0.30(7sat- 
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Figure 11: The dependence of field patterns, including the local characteristic dune density 
and size, on the sand supply rate diminshes with distance downwind: Fields appear to ap- 
proach a common attractor. A histogram of dune widths for regions 1,2,3, and 4 for the 
different sand supply rates, (A) qi, i = 0.25gsat and (B) q^i = 0.40gsat5 demonstrate that re- 
gions exhibit one, but different, characteristic size where SWq and 5Wo correspond to the 
crowded and un-crowded calving widths. As parameterized for the Tarfalya corridor, these 
results represent a 5 km by 100 km field after ~ 10000 years (~ 1800T). 
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